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Bundles of filamentous polymers are primary structural components of a broad range of cytoskele- 
tal structures, and their mechanical properties play key roles in cellular functions ranging from 
locomotion to mechanotransduction and fertilization. We give a detailed derivation of a wormlike 
bundle model as a generic description for the statics and dynamics of polymer bundles consisting of 
semiflexible polymers interconnected by crosslinking agents. The elastic degrees of freedom include 
bending as well as twist deformations of the filaments and shear deformation of the crosslinks. We 
show that a competition between the elastic properties of the filaments and those of the crosslinks 
leads to renormalized effective bend and twist rigidities that become mode-number dependent. The 
strength and character of this dependence is found to vary with bundle architecture, such as the 
arrangement of filaments in the cross section and pretwist. We discuss two paradigmatic cases of 
bundle architecture, a uniform arrangement of filaments as found in F-actin bundles and a shell-like 
architecture as characteristic for microtubules. Each architecture is found to have its own universal 
ratio of maximal to minimal bending rigidity, independent of the specific type of crosslink induced 
filament coupling; our predictions are in reasonable agreement with available experimental data for 
microtubules. Moreover, we analyze the predictions of the wormlike bundle model for experimental 
observables such as the tangent-tangent correlation function and dynamic response and correlation 
functions. Finally, we analyze the effect of pretwist (helicity) on the mechanical properties of bun- 
dles. We predict that microtubules with different number of protofilaments should have distinct 
variations in their effective bending rigidity. 

PACS numbers: 87.16.Ka,87.15.La,83.10.-y 



I. INTRODUCTION 

Bundles of filamentous polymers like F-actin form pri- 
mary structural components of a broad range of cy- 
toskclctal structures including stcreocilia, filopodia. mi- 
crovilli, cytoskeletal stress fibers, or the sperm acrosome. 
Actin-binding proteins allow the cell to tailor the dimen- 
sions and mechanical properties of the bundles to suit 
specific biological functions. In particular, the mechani- 
cal properties of these bundles play key roles in cellular 
functions ranging from locomotion 0, 0] to mechan- 
otransduction [4| and fertilization [5|. In view of this 
ubiquity, a detailed understanding of bundle mechan- 
ics is fundamental to gaining a mechanistic understand- 
ing of cellular function Quantifying the governing 
mechanical principles of these fundamental cytoskeletal 
constituents could also prove valuable in the design of 
biomimetic nanomaterials. 

In-vitro experiments recently investigated the role of 
actin-binding proteins like fascin, a-actinin and /-plastin 
in mediating bundle mechanical properties Already 
an inspection of bundle conformations from fluorescence 
microscope images makes it evident that the properties of 
the various crosslinking proteins must be quite distinct. 
While bundles formed by fascin show a compact form and 
remain straight over several microns, bundles forme d by 
a-actinin or filamin are wiggly and very lose [1, |^, ITo| . 
The mechanical properties of actin bundles formed by 



different crosslinking proteins was quantified by a fluctu- 
ation analysis 01 , which measures the magnitude of their 
thermal fluctuations. It was found that the apparent 
bundle bending stiffness can be varied over a substantial 
range by changing the type and relative concentration of 
the crosslinker. 




FIG. 1: Wormlike bundle model. We consider bundles that 
consist of regular arrangements of filament. These are as- 
sumed to be locked in place by crosslinking proteins. When 
the bundle bends and twists in space, the filaments start to 
slide along each other. This effect leads to shear deformation 
in the crosslinks. 
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These intriguing mechanical properties can be under- 
stood in terms of the wormUke bundle model (WLB), 
which describes bundles as an assembly of semiflexible 
filaments interconnected by crosslinking proteins pj].[T^: 
for an illustration of the bundle architecture see Fig. [11 
Unlike the standard wormlike chain model (WLC) [IJ, 
[T^ . the wormlike bundle model (WLB) exhibits a state- 
dependent bending stiffness [11| that derives from a 
generic competition between the bending and twist stiff- 
ness of individual filaments and their relative motion me- 
diated by the stiffness of the crosslinkers. An important 
aspect of the WLB model is that crosslinks may be very 
efficient in constraining the lateral excursions of filaments 
within the bundle but much less so in inhibiting their ax- 
ial motion. This is possible as the relative axial sliding 
of two crosslinkcd filaments probes not only the clastic 
properties of the crosslinking protein but also those of 
the binding domain at which the protein is attached to 
the filament. The latter may be of quite different rigidity 
than the protein itself. Unfortunately there are no single- 
molecule experiments yet which would quantify the me- 
chanical and binding properties of actin crosslinking pro- 
teins attached to a pair of F-actin filaments [G^I ■ In the 
WLB model the mechanical properties of the crosslinks 
are described by a single shear-stiflFness fcx for the relative 
sliding of the constituent filaments. 

An important finding within the wormlike bundle 
model is that the mechanical properties of bundles can 
be classified into three distinct bending regimes that are 
mediated by both crosslink type and equally importa ntly 
by bundle dimensions, namely diameter and length 
Taking into account the mechanical properties of fila- 
ments and crosslinkcr at a microscopic level is a virtue of 
the model making it quite verstilc. It reduces to a well de- 
fined continuum limit but is equally applicable to bundles 
with as few as two filaments. This microscopic perspec- 
tive may prove a valuable starting point to address more 
complex questions related to bundle mechanical proper- 
ties. This may include problems such as disorder, lattice 
defects [l^, or filament fracture [T6| . 

While previously we have provided a formulation only 
for plane-bending in two dimensions here we give 
a full derivation of the WLB Hamiltonian in three spa- 
tial dimensions. This includes bending as well as twist 
deformations. Wc explore the predictions of the WLB 
model for experimental observables such as the tangent- 
tangent correlation function and dynamic response func- 
tions. Moreover, wc discuss the effects of different bun- 
dle geometries on their mechanical properties. In this 
respect, we will view microtubules as a bundle of (proto- 
)filamcnts arranged on the surface of a cylinder; compare 
Fig. [21 This shell-like bundle architecture is contrasted 
with a uniform distribution of filaments as found for F- 
action bundles EBl ■ Finally, we discuss how hclicity 
influences bundle mechanics. 




FIG. 2: llustration of the bundle geometries considered. The 
position of the kth filament in the cross-sectional plane {{y, z)- 
plane) is given by R^. F-actin bundle architecture (left): fil- 
aments are arranged on a square lattice. Microtubule archi- 
tecture (right): filaments are arranged on the surface of a 
cylinder. 



II. MODEL DEFINITION 

We consider bundles of length L that consist of N par- 
allel filaments. While the filaments may form a disor- 
dered lattice structure, we will focus our attention here 
to the cases of regular arrangements. In particular wc 
will treat in detail the square-lattice and the cylindrical 
tube (see Fig. [2]). 

Each filament is modeled mechanically as an exten- 
sible worm-like polymer with stretching stiffness fcs, 
bending rigidity and twist rigidity kj. Filaments 
are irreversibly crosslinkcd to their nearest-neighbors by 
crosslinks with mean axial spacing 5. The crosslinks 
are modeled to be compliant in shear along the bundle 
axis with finite shear stiffness , and to be inextensible 
transverse to the bundle axis, thus constraining the in- 
tcrfilament distance, 6, to be constant (see Fig. [2). This 
assumption, which neglects crosslink stretching deforma- 
tions, is based on the recognition that the shearing mode 
involves deformation of the crosslink and its binding do- 
main to the filament. The resulting stiffness may indeed 
be much lower than that of a crosslink in isolation. 

Filament stretching is characterized by the axial dis- 
placement Ufe(s) of filament k at axial position s along 
the backbone. To describe bundle bending and twisting 
wc define {di,d2,d3} to be a material frame fixed to 
the bundle central line at each arclength position s. The 
vector d3 = t is the tangent to the space curve traced 
out by the central line, while the two vectors di and d2 
lie within the cross-section of the bundle. The position 
of each filament in the cross-section is parametrized in 
terms of a vector Rfc(s) = Afcdi(s) -|- i?fed2(s), where 
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Ak and Bk are the material-frame coordinates of fila- 
ment k; they are constants independent of arclength s 
and deformation of the central line (see Fig. [2|). As one 
moves along the bundle backbone the material frame ro- 
tates according to Frenet-Seret equations, dgd = x d. 
The rate of rotation is given by the generalized curva- 
tures n — (r^i (s), r22(s), ri3(s)), which, in addition to 
the axial displacement u^, represent the basic kinematic 
degrees of freedom of the bundle. 



A. The WLB Hamiltonian 

Neglecting all nonlinear effects we arc now going to 
develop a simple expression for the bundle energy which 
is harmonic in its degrees of freedom, axial displacement 
Uk and generalized curvatures Qa. 

This WLB Hamiltonian consists of three contributions, 
ffwLB = Ho + ^^strotch + i?shoar- Thc first tcmi corrc- 
sponds to the standard wormlike-chain Hamiltonian 

N 

Hq = ^ ds (fi? + nl) + Ktn"^ . (1) 

^ Jo 

Writing this we assume that each filament follows effec- 
tively the same space-curve as the center-line. While, 
in general, one should account for the curvatures, fifc, 
of each individual filament separately, this would only 
lead to correction factors that can be neglected for our 
purposes. Consider, for example, planar bending of the 
central line, fii = l/p, where p is the radius of curva- 
ture. The filaments that lie at a distance R away from 
the central line naturally have a different radius of cur- 
vature, p ± R, and thus a different bending energy. The 
magnitude of the correction term relative to fi^ is small 
since it scales as i.e. we assume the typical cur- 

vatures to be smaller than the bundle radius (for more 
details see Appendix |B|) . 

As to the twist degree of freedom, in Eq. ([T|) we do 
not allow for the possibility of relative twisting of the 
individual filaments (see Section |V] and Appendix [B] for 
a discussion of this effect). We assume that twist is only 
due to thc "bundle-twist" il^ of the central line. This 
assumption is reasonable in tightly bound bundles, where 
the filaments are connected by many crosslinks. In this 
state the filaments and their relative orientation can be 
assumed to be locked-in by the crosslinker. 

The second term in the Hamiltonian, iJstrctch = 
Sfc -^stretch' accounts for filament stretching. It depends 
on the difference in axial displacement, Uk, between two 
crosslinks at arclength positions Si and Si + 5, respec- 
tively. 



where we have performed the continuum limit — > 
J ds/S to arrive at the second line. The spring constant 
ks (S) is the single filament stretching stiffness on the scale 
of the crosslink spacing S. 

The particular form for ks depends on the system 
under consideration. For high crosslink concentrations 
(small S), the segment behaves as a homogeneous elas- 
tic beam, characterized by a Youngs modulus E and 
^boam ^ Eb'^/d. The combination ksS that enters the 
Hamiltonian is independent of 6, as it should: the me- 
chanical stretching stiffness of a beam cannot depend on 
the properties of the crosslinks. 

For small concentrations of crosslinks (large 5) entropic 
effects become relevant and the stretching stiffness is that 
of a thermally fluctuating wormlike chain with persis- 
tence length Ip. In this case one has fc™*'' ~ Kt,lp/6'^, 
which implies that thc combination kgS ^ does de- 
pend on thc crosslink spacing S. This is related to the 
fact that the formation of a crosslink suppresses ther- 
mal undulations (reduces entropy) and thus increases the 
entropic stretching stiffness. Equating both stretching 
stiffnesses, fc^"™ ~ /c™'"' one finds the critical crosslink 
concentration, ~ b'^lp-, at which the cross-over from 
enthalpic to entropic elasticity takes place. 

In the case of microtubules, which will be treated in 
Section lTVl the crosslink spacing S is given by thc tubulin- 
sizc and the stretching stiffness is modeled as for an elas- 
tic beam. 

The third energy contribution, i?shoar = Life -^shcar' 
results from thc crosslink-induced coupling of neighbor- 
ing filaments. The relative axial motion of a filament 
pair {I, k) at a given point of the backbone is described 
by the crosslink shear displacement, which is the sum of 
a geometric contribution, A/^ , and the relative stretching 
of neighboring filaments, Auik = ui — Uk- The geomet- 
ric part results from the arclength mismatch between the 
two filaments, induced by a bending and twisting of the 
bundle central line. As in Eq. ^ we first write thc shear 
energy as a sum over all crosslink positions Si and then 
perform the continuum limit, to get 

HiU = ^^[Aa-(s.) + AK,,(..)]' 




L 

ds {Aik + Auikf , (3) 



where kx is the shear stiffness of the individual crosslink. 

For any bundle deformation, the associated value of 
A/fe can be compensated for by stretching the filaments, 
making the shear energy vanish when Auik — — A/^. At 
the same time, however, this would increase the stretch- 
ing energy, which may be unfavorable if thc stretching 
stiffness kg is rather large. For deformations on the scale 
of the bundle length {u' ^ u/L) the ratio of both en- 
ergies gives the important parameter a = ky,l? jksS^ ^ 
which quantifies the relative strength of both deforma- 
tion modes [ll| . 

As a final ingredient to the model we need to calculate 
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the dependence of A^- on the bundle curvatures, fta- 
Without going into the details of an explicit derivation, 
we here just give the resulting expression. For more de- 
tails we refer the reader to Appendix [C] The special case 
relevant for the microtubule geometry is also dealt with 
in some detail in [l9l |. To linear order in Qa we find 



bik cos aiklyik^S" dtfliit)) (4) 



^Ik 



-bik sin aik (^ik^^ - j dt^i{t)^ , 

where we defined bik — |b;fc| as the distance between the 
filament pair and aik as the angle of hik with respect 
to the z-axis. Furthermore, (yikTZik) are defined as the 
cross-sectional coordinates of the midpoint between the 
filament pair. 

In contrast to Hq, which is an expansion in the general- 
ized curvatures Qa-, the shear energy iJshcar is a function 
of the integrated curvatures, since A ~ & f2 ^ bLfl. For 
terms beyond the harmonic contribution to be negligible, 
one has to assume that the shear displacement is suffi- 
ciently small, A ^ a, where a is some microscopic length- 
scale related to the size of the crosslink. In terms of the 
bundle curvatures this implies Lfl <^ a/ R ~ 1, which is 
much more restrictive than the range RVl <^ 1, over which 
Hq can be approximated by a harmonic form. As a con- 
sequence the bundle is only allowed to make small excur- 
sions from its initial state, an assumption which is usu- 
ally well satisfied in bundles of stiff polymers like actin or 
microtubules, but certainly breaks down under extreme 
loading conditions (e.g. to describe post-buckling) or in 
more fiexible objects like DNA. 

With this "weakly-bending" assumption we reformu- 
late the generalized curvatures in terms of the lab-frame 
Euler angles 0, ^ 



ill 

^3 



dcj} . , . ^ de 

— — sm ip sm tl + — cos ip 

OS OS 

, . n 99 
— — cos ip sm a — — sm %p 
OS OS 

dtjj d(j) 



ds 



ds 



cosf 



As reference state we take (po = tt/2, 9q = 7r/2 and 
ipo = sujQ, which corresponds to a straight, but pre- 
twisted bundle that points along the cc-axis. For small 
excursions around this reference state we can linearize 
the equations such that 



ill = -^sn\ipQ + —cosTpQ 
as OS 



(5) 



"2 = t;- cos ■(/'o - TT ■(/'o 
as OS 

as 

and the angles are now measured relative to the reference 
state. As expected the prctwist i/iq leads to a coupling 



of the angles 6 and (j). In the following we arc primarily 
concerned with the case of vanishing pretwist. Then, the 
coupling terms vanish and we can simply set 

We will come back to the case of pretwist in Section |Vl 

B. Examples for the arclength mismatch A 

For the purpose of illustration we provide some ex- 
amples of how the shear displacement A depends on the 
geometry of the bundle cross-section and the deformation 
of its central line. 

If the bundle consists of only two filaments [l^l (ge- 
ometry of a ribbon) we have y = z = a — tt/2 as the 
central line and the mid-line between the two filaments 
are identical. In this case A simplifies to 



A 



(6) 



which is illustrated in Fig. [3] Note, that twist (f^a) does 
not contribute, as both filaments twist around the central 
line symmetrically. 




FIG. 3: Illustration of crosslink shear deformation for the case 
of a two- filament bundle (filaments in red). Bundle deflection 
through the angle 9 leads to the arclength mismatch, A = b9. 
The filaments have to stretch the relative amount u\—U2 = h9, 
in order to keep the crosslink (dashed line) undeformed with 
zero shear energy. 

A particular case of this two-filament bundle has been 
considered in a set of articles [23l [2^ . [2^, where the 
crosslinks are assumed to be rigid with respect to shear, 
i.e. fcx oo. To satisfy a vanishing arclength mismatch 
A one thus requires r2i = — d2-t = 0. This means that 
the unit vector d2, which points from one filament to 
the other, must not rotate in the direction of the tan- 
gent t. In other words, d2 must equal the the bi-normal. 
Thus the ribbon orientation is completely specified by 
the space curve traced by the central line. 

A second example is the axoneme in eukaryotic flag- 
ellae [l^. There, filaments (microtubule doublets) are 
arranged on the surface of a cylinder, just as the protofil- 
aments in a single microtubule. Switching to polar coor- 
dinates in the cross-sectional plane (see Fig. (y, z) 
(R, if), a = (fi + tt/2, we find 



A 



-bRip' + b(t) sin (f + b6 cos ip . 



(7) 
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action of a bending moment. In the present formula- 
tion of the model these nonlinearities are not accounted 
for. For a discussion of potential modifications to include 
cross-sectional deformations, we refer the reader to the 
/ outlook section at the end of this article. 




FIG. 4: Illustration of crosslink shear deformation for a 
twisted microtubule. Indicated (in red) is a filament pair that 
winds around the microtubule cylinder (radius R) taking an 
angle R^ls with the cylinder axis t. The resulting arclength 
mismatch is given by A = bR^ls. 



A similar expression, disregarding the possibility of twist- 
ing, has been given by Mohrbach et al. [26|]. The struc- 
ture of the second and third term is the same as in 
Eq. ([6]), additionally taking into account the modified 
orientation in the cross-sectional plane, as described by 
the angle ip. The origin of the first term is illustrated in 
Fig. m and elaborated on in Appendix [C] 

Consider now the case of planar bending. This will 
make the connection to continuum elasticity particularly 
clear. Under planar deformation the bundle is described 
by the one variable 6*' = such that the shear defor- 
mation is simply A = hO{s) ~ hdxUy. Here, Uy is the 
displacement of the bundle transverse to the bundle axis 
(in the y-direction) . This latter form makes clear that 
the shear deformation represents one part of the strain 
tensor component e^y = \{dxUy -\- dyUx). The second 
part, dyUx ^ bAuik is the continuum version of relative 
filament stretching (see Eq. 

The elastic symmetries relevant to the WLB model 
depend on the arrangement of the filaments in the cross- 
section. If there is rotational symmetry with respect 
to the bundle axis, one speaks of transversly isotropic 
elastic bodies [l^]. While, in general, this has five first 
order clastic constants, our model has only two, aug- 
mented by the (second order) bending/twisting elasticity 
of the individual filaments which is not accounted for in 
continuum elasticity. The simplification arises from as- 
suming transverse inextensibility as well neglecting cross- 
sectional shape changes. The latter is, for example, im- 
portant in the failure of hollow tubes under bending. One 
relevant effect is the Brazier effect (2^, which describes 
the increasing ovalization of the cross-section under the 



C. Definition of effective bending and twist 
rigidities 

It should be clear from the way the Hamiltonian was 
derived that the model is applicable to bundles with ar- 
bitrary (ordered/disordered) arrangements of filaments 
in the cross-section. In the remainder of this article we 
will focus our attention to bundles with highly symmetric 
cross-sections, where the filaments cither form a rectan- 
gular array or a hollow tube. 

In view of recent experiments probing the mechanical 
or statistical properties of individual bundles in vitro 0] , 
we head at a description of the bundle in terms of effec- 
tive bending and twist rigidities. These are defined with 
respect to the standard worm-like chain model. To ar- 
rive at the proper expressions we have to integrate out 
the internal stretching variable u, which in general is not 
observable in experiment. 

To show how this works, we symbolically write the 
partition function a,s Z = J D(f)Z{(p) where (j) sig- 
nifies the set of Euler angles 4>{s),9{s),tp{s). The 
constrained partition function then reads Z{(p) = 
/ D{u} exp(-/3i7({0, u})) = exp{- pW {(f))) . 

The integration over the w-variables can easily be per- 
formed. As the Hamiltonian is harmonic we are left with 
only Gaussian integrals, which are evaluated in Fourier 
space. The resulting potential of mean force W{4>) can 
be written in the form of a wormlike chain Hamiltonian 

Wi<t>n) = ^Y.1n [^BinMl + el) + KT{n)^l^l] (8) 

n 

with effective bending and twist rigidities KB(n) and 
KT{n), respectively. We note that in the symmetric sit- 
uations considered here there is no coupling between the 
different deformation modes bending and twisting (see 
Appendix [A]) . In contrast to the usual WLC, the effec- 
tive bend and twist rigidities are in general dependent on 
the mode-number n and thus on the wavelength of the 
deformation. This effect and the discussion of its conse- 
quences is the central topic of the remaining sections. 

III. F-ACTIN BUNDLE ARCHITECTURE 

In the following sections we will focus our attention 
to bundles with N = (2M)^filaments that form a rect- 
angular array (see Fig. The angle a that specifies 
the orientation of the filament-pair in the cross-section is 
then a = 0, 7r/2 as filaments are either arranged along the 
y- or the z-axis. As mentioned above the different defor- 
mation modes decouple in harmonic order. We can thus 
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investigate bending independently from twisting. Also 
the two space-directions decouple and we can reduce the 
model to an effective two-dimensional description [l^ . 



A. Effective bending rigidity 



The shear Hamiltonian reduces to 



M-1 

Jo ,r,n 



Ufc + bO)^ , (9) 



where we used Eq. (|4]) with (p = ip = 0. By following 
the recipe outlined above we eliminate the axial strain 
variable Uk- By approximating Uk by a linearly increas- 
ing function of k (see discussion below) we arrive at the 
following result for the effective bending rigidity [H^l as 
defined in Eq. © 



KB{n) = NKb 



12kb 
N~ 1 



Here, we have defined a characteristic wavelength 



2M 



{2M - 1) 



(10) 



(11) 



and a dimensionless bending rigidity ki, — Kf,/(fcs(56^). 
In terms of the quantities A and kt the previously de- 
fined a = kxL^/ks6^ can be rewritten as a ~ kt,{L/X)'^. 
If the filaments behave as homogeneous elastic beams, 
kb is just a number independent of bundle geometry or 
crosslink spacing. For any numeric computation we will, 
for specifity, assume that kb ~ 1/12, which corrrcsponds 
to beams with square cross-sections 0] . 

The characteristic feature of Eq. pUj) . is the wave- 
length dependence (see Fig. [5]). For wavelengths in 
the interval ^ qnX <C 1 the bending stiffness de- 

creases as KB{n) ^ fcx^n^- 111 -fl^f- [ill have termed 
this the intermediate or shear dominated regime as the 
bending rigidity is proportional to the shear stiffness of 
the crosslinks. It is in this parameter regime that the 
bundle behaves qualitatively different than either a ho- 
mogeneous beam (obtained in the " fully coupled" limit of 
q„A <C l/x/iV) or an assembly of "decoupled" filaments 

(Va> 1). 

Another important feature of Eq. p^ . which is in- 
dependent of the specific (/-dependence, is the ratio of 
maximal to minimal bending rigidity, r = Kmax/'Smin = 
1 -I- {N — l)/12kb. This only depends on the number of 
filaments and the dimensionless bending rigidity kb- 

While integrating out the stretching variables Uk can 
be performed exactly, Eq. (fTO|) is based on the addi- 
tional assumption that axial strains are linearly increas- 
ing through the cross-section, Uk ~ Au ■ {k + 1/2). The 
exact profile for is calculated in the appendix and dis- 
played in Fig. ([S]); compared to the linear profile it shows 
an enhancement of strain towards the bundle periphery. 




FIG. 5: Effective bending rigidity, Eq. (|10p . as a function 
of mode-number, qX, and for the set of bundle sizes A*" = 
4, 9, 16, 25 (from bottom to top). In the fully-coupled and the 
decoupled regimes (corresponding to small and large q) the 
bending rigidity is constant. At intermediate values of q the 
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FIG. 6: Dependence of axial strain on k, the distance 
from the bundle central line. Decreasing the dimensionless 
shear-stiffness a = k^L/^ /ksS^ the strain is reduced but not 
in a linear fashion. The outer layers of the bundle remain 
stretched stronger than the inner ones. 



However, the ensuing value for the bending stiffness is 
largely insensitive to the linear approximation [g^I- We 
speculate that the nonlinearities in the axial strain may 
eventually be important for nonlinear material proper- 
ties, as for example strain induced rupture. The in- 
creased strain in the outermost filaments brings them 
closer to their threshold for rupture and thus makes them 
more susceptible to this mode of failure. 

In order to make contact with continuum models for 
beam bending we perform a continuum limit, by letting 
^ oo but keeping the bundle aspect ratio D/L ^ 
hM / L constant. Then, bundle length L has to grow with 
M as L{M) ~ M . In particular, this implies that fewer 
and fewer modes n belong to the decoupled regime (where 
q^^ <C A). Eventually, this regime, where filaments bend 
independently (/{^(n) w Nub) is no longer accessible. 
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In effect this means that the bending stiffness Kb of the 
individual filaments can be neglected, just as in "nor- 
mal" continuum elasticity, where higher order gradients 
{0{6')) are not accounted for from the start. 

In this continuum limit the result from the lineariza- 
tion assumption, Eq. pO|), reduces to the Timoshenko 
model for beam bending [30[, which was recently used 
to inter pre t bending stiffness measurements on micro- 
tubules |3ll . [3^ and carbon nanotubc bundles [s^ , 



TIMO 



(n) 



1 + {qnDYE/l2G 



(12) 



To allow for comparison with continuum elasticity wc 
have used the expressions kgS = Eb^ and ki, — Eb'^/U 
applicable for homogeneous beams of square cross-section 
and defined the shear-modulus G = k^/S. 

On the other hand, one can equally derive a continuum 
limit from the exact expression for kb(7i) (as presented 
in Appendix I A 1|) . This gives 



,CONT 



in) 



1 - 



tanh(g„i^/E74G) \ 
q.nD^/Wj4G J 



{qnD)^E/l2G 

(13) 

Both expressions, Eq. p2)) and Eq. (|13p arc compared 
in Fig. ([7]). One infers that the approximation (Tim- 
oshenko theory) overestimates the exact solution by no 
more than 6%. The difference can partly be compensated 
for by introducing a shear-correction factor (« 1.2) in the 
denominator of Eq. (fT2|l . 




FIG. 7: Continuum limit of effective bending rigidities as a 
function of a: = (qnD)^ E/12G. Comparison of exact solution, 
Eq. (|13p . with approximation, Eq. (|12p . Inset: Relative error 
between both expressions. The approximation over-estimates 
the bending stiffness by maximal ~ 6%. 



B. Tangent-tangent correlation function 

In this section the implications of a mode-number de- 
pendent bending rigidity is further elaborated by dis- 
cussing the concept of the persistence length. The per- 
sistence length of a single WLC may be defined in terms 



of the competition of bending and thermal energies, 
Ip = n/kBT. With this definition, bending rigidity and 
persistence length are basically identical. In the frame- 
work of the WLB this would lead to a mode-number de- 
pendent persistence length lp{n) = KB{n)/kBT. 

For the WLC the persistence length is also the length- 
scale over which the tangent-tangent correlation function 
decays 



(t(s)t(O)) =cxp[-s/y 



(14) 



This simple exponential form is no longer valid for the 
WLB as can be illustrated by considering planar bending 
with (/) = 0. Then the tangent-tangent correlation func- 
tion is easily inferred from the angular fiuctuations 



(t(s)t(O)) = exp 



with 



1 



A-l 



A 



1 



(15) 



s^/A/X 



AVA 

(16) 

and A = 1 + 12ki,/{N — 1). Forcing such a complex 
expression into the form given by Eq. (|14p implies an 
arclength-dependent persistence length, lp{s). At short 
distances, one recovers the decoupled regime and lp{s) = 
Nip, while at long distances, lp{s) = N'^lp as found in 
the fully coupled regime. 

Note, that there is no immediate relation between this 
lp{s) and the lp{n) defined above. The Fourier-transform 
of lp{n) is, instead, given by the following expression 




(17) 



This quantity appears in the elastic energy expressed in 
real-space as 



i^WLB = kBT J 9'{si)lp{si - 82)9' {s2)dsids2 , (18) 

which is a non-local function of arclength. The length- 
dependence obtained here is markedly different from that 
found in the correlation function, Eq. ([T6|). While lp{s) 
is constant at large distances, lp{s) decays exponen- 
tially and vanishes over the length-scale \/\/ A — 1 ^ 
Xy'N/kb, which corresponds to the onset of the fully- 
coupled regime. 

A similar nonlocal energy function is obtained when 
considering the fluctuation properties of elastic mem- 
branes. In-plane shear- and compression-modes lead to a 
renormalizcd bending rigidity for the out-of-planc fluctu- 
ations [HjII^. In contrast to bundles, however, there the 
kernel is long-ranged lp{s) ~ 1/s, which asymptotically 
leads to a flat membrane phase. 

Concluding this section, we find that it is impossible 
to speak of a single persistence length without specifying 
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the precise experimental conditions as well as the observ- 
able under consideration (here tangent-tangent correla- 
tion function). The WLB model presents a framework 
within which a length-dependent bending stiffness can be 
understood. However, we would like to emphasize that 
the fundamental quantity is the (jn-dependent KB^n) as 
presented in Eq. (fTUl) or below in Eq. It depends on 
wavenumber in a universal way independent of the type 
of measurement. The dependence on bundle length, in 
contrast, arises through a specific transformation to real 
space and may produce different results depending on the 
observable under consideration. 



Frequency-dependent correlation and response 
functions 



In our previous publications [13] we have calculated 
several thermodynamic observables and showed that the 
mode-dependent bending stiffness of a WLB may lead 
to drastic modifications of their scaling behavior. Here, 
we widen the scope of this analyis by discussing dynamic 
observables. In analogy to the usual overdamped dynam- 
ics of a WLC we can discuss the dynamics of a WLB by 
substituting the mode-dependent bending stiffness in the 
standard Langevin equation of motion for the transverse 
bending modes r±(qmt) 



' dt 



(19) 



With this one obtains for the reduced correlation func- 
tion 



at) 



L- 



1 - e"*/^" 



(20) 



with the relaxation times t„ = C,/ KB{n)q^. For WLCs 
(constant kb = k) one finds a scalin g regirn e at times 
t < ri - (L^/k, where C{t) - t^/^ [31^33, |3|. 

For WLBs, on the other hand, one has to use the q- 
dependent effective bending stiffness, Eq. pU]) . instead. 
As the term {qX)~^ multiplies the q* contribution in 
Eq. (PO)) one finds the correlations to grow in time as 
CwLB(i) ^ t^^^- Of course, this result is only vahd as 
long as the gA-term dominates the effective bending stiff- 
ness, i.e. as long as the bundle is in the intermediate 
regime. In fact, there will be a complex cross-over see- 



C{t) 



(Nip)-' 



TV 

Kb 

N 



c 

1/4 



1/2 



NKbt 



3/4 



t<ii 

tl < t < t2 
t2 < i < ^3 

(21) 

Kb 



N 



with the cross-over times ti = (X /NKb and ^2 = ^i 

governing the cross-over from the decoupled to the inter- 
mediate and the fully-coupled regimes. At times larger 

Kb 

than t-i = t-\ — r the correlation function saturates. 

N X^ 

From the correlation function one can furthermore cal- 
culate the response function X-L measuring the linear 
response to transverse forces. Using the fluctuation- 
dissipation theorem and the Kramers-Kronig relations 
one finds 



1 



L K.B{n)qf^ - iC,uj 



(22) 



This contrasts with the response function X|| for stretch- 
ing forces 



_1 ^ ^ (23) 

4^ lp{n) KB{n)qi - iQuj/2 ' 



which is sometimes taken as a starting point to determine 
the high-frequency response of entangled solutions of stiff 
polymers ^39, 40, 41, it2, 4^], The transverse response, on 
the other hand, has been argued to relate to a "microrhe- 
ological modulus" 0, that is measured locally by 
imbedding probe beads into the network. Unfortunately, 
for a WLC both functions are hardly indistinguishable, 
and only differ by the constant factor X||/X-L = L/lp. 
The high-frequency behavior in both cases is x ~ uj^^^^. 

In the case of a WLB, things are somewhat different, as 
the additional factor L/lp{q) in the longitudinal response 
function not only changes the prefactor but also modifies 
the functional form with respect to frequency (see Fig.[H]). 

For similar reasons as in the discussion of the corre- 
lation function one expects an intermediate regime with 
ll!^^/^, at least in the transverse response. This should 
lead to measurable signatures in microrheological exper- 
iments on bundled F-actin systems. Due to the addi- 
tional q-dependence in the denominator, the longitudi- 
nal response, X|| shows a smooth cross-over between the 
two asymptotic regimes of fully-coupled and decoupled 
bending, as explained in the figure. 
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FIG. 8; Comparison of transverse response function x±{^) 
(top) with longitudinal response function X\\{^) (bottom) for 
a bundle of A'' = 20 filaments (in units of Nki,/L^ and as- 
suming LksT /N Kb = 1). Frequencies are plotted in units of 
(L^/Nnb. The different curves correspond to different val- 
ues of X/L. The two asymptotic scaling regimes with uj~^'''^ 
correspond to the decoupled (top) and fully-coupled regime 
(bottom), respectively. The intermediate is sharper in 

X± than in x\\ even though the latter function shows overall 
a stronger variation with frequency. 



D. Effective twist rigidity 



NKt 



N - 



1 /At 
b 



(25) 



Here we have defined A* := V2W(2AV^T) V^tV^^, 
which is similar to A defined above, with Kf, substituted 
by K(. Unlike A, however, it multiplies the constant 
length b, the lateral distance between neighboring fila- 
ments. As anticipated, there is no mode-number depen- 
dence. 

By introducing bundle diameter D and effective shear 
modulus G (as in the discussion leading to Eq. (H^)) the 
second term takes a form well known from continuum 
theory, - GD* Thus, we find an effective twist 

rigidity, Eq. that just describes the simple additive 
superposition of two contributions, shear-induced rigid- 
ity (second term) and twist rigidity of the individual fil- 
aments (first term). In the continuum limit N —^ oo 
the latter contribution can naturally be neglected as it 
only grows with TV as compared to the N"^ increase of 
the shear-induced rigidity. 

Now we allow for finite axial displacements Ufc. This is 
commonly referred to as cross-sectional warping; under 
twist deformations the bundle cross-sections do not stay 
plane but deform and acquire a curvature. Just as in the 
case of the bending rigidity, the exact solution for the 
twist rigidity only differs little from the foregoing sim- 
plified analysis. Some details about the derivation are 
presented in Appendix I A 21 The resulting axial displace- 
ments Uk can be found in Fig. The classic solution 
of Saint- Venant [23| is approached closer and closer for 
increasing the shear-stiffness fcx . 



We now turn to the discussion of the twist mode. In 
this case the Euler angles (p and 6 in Eq. ^ are zero such 
that the shear Hamiltonian reduces to 

ifshear = ^ / ^ ( Au,fc + bdlk^'f , (24) 
■^0 Ik 

where we defined the geometric factor dik = yik cos aik — 
zik sin aik ■ 

As may already be apparent from comparing Eq. ([9|) 
with Eq. the stretching deformation u couples differ- 
ently to twist {k-^uij}') as to bending {kxu6). The differ- 
ence being the additional derivative occuring in Eq. ([24|) . 
Effectively this means that the resulting twist rigidity will 
not depend on mode-number g„ but receive a constant 
correction to the single filament value Kt- 

Let us first assume that the filaments are inextensible. 
Then the w-terms in the Hamiltonian identically vanish 
and the effective twist rigidity can simply be read off from 
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FIG. 9: Cross-sectional warping Uk {k — —M + 1 . . . A4 , cor- 
responding to one row of the rectangular array) for differ- 
ent non-dimensional crosslink shear stiffness a = kxL'^ /ksS^ . 
The classic solution for a beam (dashed line) is approached as 
fcx oo. In the opposite limite, fcx 0, there is no warping 
and u = 0, as filaments remain uncoupled in this limit. 
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IV. MICROTUBULE ARCHITECTURE 

In this section wc want to turn our attention to the 
case of microtubules, which we model as bundles with 
filaments arranged on the surface of a cylinder. For the 
time being we assume that in the groundstate the micro- 
tubule is untwisted such that the N protofilaments are 
oriented along the cylinder axis. This assumption is in 
fact only valid for microtubules with N = 13 protofila- 
ments [i^ll^]. This class, nevertheless, seems to be the 
most frequent in in-vitro polymerization experiments. 

There is an ongoing debate in the literature about the 

tendence of microtubule bending rigidity on length [bH , 
m, m, M, \M,\M,\EM- Early buckling experiments 
indicated a length-dependence [48j . while an improved 
version of the same experiments later gave a negative 
result [4^ . Brangwynne et al. [50i | performed a mode 
analysis of microtubule contours. Their data is compati- 
ble with a length-dependent bending rigidity but the au- 
thors vote for a cautious interpretation of the results in 
view of the close proximity to the noise level. Recently, 
experiments by Pampaloni et al. (3^ measured the trans- 
verse fluctuations of grafted microtubules to establish an 
increasing persistence-length for microtubule lengths up 
to L ^ 23/im. Using a high-precision tracer technique, 
Taute et al. [IH analyzed shorter microtubules and found 
the persistence length to level off at Z™'" « 580//to for 
lengths shorter than L sa 5/im. Similar values have been 
obtained in Ref. ^ (Z^f'" = 90/^m) and in Ref. H 
(?™'" = 240^m), and explained with the help of the WLB 
model. On short length-scales (in the decoupled regime) 
the effective bending rigidity is constant because it re- 
flects the stiffness of the individual proto-filamcnts. In 
contrast, Kis et al. (sij have found a decreasing bending 
stiffness even for microtubule sections as short as several 
hundred nanometers. 

While the discussion about these partially conflicting 
measurements is ongoing, we would like to point out that 
different techniques do not necessarily have to come to 
the same conclusion, once the idea of the bending rigid- 
ity as a fundamental material parameter is given up. In 
the context of the WLB model the fundamental quantity 
is a mode-dependent bending rigidity. As discussed in 
Section llll Bl any dependence on bundle length is a "sec- 
ondary" effect that will depend on the observable under 
consideration. 

By using the same procedure as in the case of the 
rectangular bundle we find for the microtubular bend- 
ing rigidity 



KBin) = NKb 



8kb 



(tt/N) 



, (26) 



w here the r elevant length-scale A is now defined as A = 
y^KbS/ky^b^ . Note, that this expression results in the 
bending stiffness in the fully-coupled regime (where the 
microtubule behaves as a simple beam) to scale as k ^ 



N^, in contrast to bundles with a homogeneous cross- 
section (as the one discussed before), where k ~ N"^. 

With Eq. (|26p direct comparison with experimental 
data can be made. In particular, the ratio of maximal 
to minimal bending rigidity can most easily be deter- 
mined, r = 1 + sin"^(7r/7V)/8K6 = 1 + 2 / [t: / N) . For 
the latter equality wc assumed the protofilamcnt to have 
a circular cross-section (Kb = 1/16). For microtubules 
with = 13 protofilaments this results in a universal 
ratio r « 35. 

In this way, using the the range of values I™™ w 
0.1 .. . 0.6mm (which we assume to represent the protofil- 
amcnt stiffness) one can estimate the maximal persis- 
tence length of long microtubules to be in the range 
K, 3.5 . . . 21mm. Compared with experiments |32l . 
[13, UHl these values seem to be too large. Given the large 
error-bars in any of the mentioned experiments this cal- 
culation may, nevertheless, be acceptable. Furthermore, 
as we will discuss below, microtubule helicity can provide 
a mechanism to reduce the apparent persistence length 
by reducing the effect of shear-induced coupling. For 
given microtubule length the persistence length is then 
predicted to decrease with increasing helicity - thus im- 
proving the comparison with experiment. 

Finally, let's turn to the case of pure twist deforma- 
tions. Due to the symmetry of the circular cross-section 
no warping is possible. We find for the microtubule twist 
rigidity, similar to Eq. ((25)) . 



tair^{7:/N) fXt 



(27) 



where Af := y^KtSjk^l^. 



V. PRETWISTED BUNDLES 

In the previous sections we have restricted our at- 
tention to non-helical bundles and assumed that in the 
ground-state the fllaments point along the bundle axis. 
As a final application of our model, we will here discuss 
the question of helicity, or prctwist, and its influence on 
bundle mechanics. This aspect is important not only for 
some types of microtubules but also for F-actin bundles 
and is reflected in the role that helicity plays in provid- 
ing an explanation for the existence of a (thermodynam- 
ically) preferred bundle size 

The discussion of pretwisted bundles proceeds in two 
steps. We first assume the bundle to form with all 
filaments straight. Starting from this reference state 
crosslink binding may add a driving force for twisting 
the bundle if the straight state docs not allow for opti- 
mal accessibility of the crosslink binding sites. Another 
effect may be the helicity of the filaments themselves that 
favor a twisted bundle over an untwisted one. 

Without elaborating on the precise mechanism that 
leads to pretwisted bundles, we incorporate bundle 
prctwist by substituting the generalized twist-curvature 
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1 10 100 

FIG. 10: Effective bending stiffness (in units of Nki,) as a 
function of shear stiffness kx (in units of KbS/b^L^) of a bun- 
dle of four filaments under a tip force F. The bending stiff- 
ness is calculated from the determined end-deflection y{L) by 
Keff = FL'^ /12y{L). The full curves are for finite pretwist 
looL/tt = (0.25, 0.5, 1, 2, 3), while the dashed curves represent 
the limiting cases of zero and infinite pretwist, repsectively. 



by ri3— cjo- In doing so the new energetic ground-state 
is at 17i = = and Vl^^ — ujq . To obtain the effective 
bending rigidity we linearize the curvatures around this 
ground-state, as performed in Eq. 

Inserting this result into the shear deformation, 
Eq. ([D), one finds terms like h J cos{iPq) (f)' ds which de- 
pend nonlinearly on arclength s. Thus a transformation 
to Fourier-space is no longer helpful as different modes 
would remain coupled. We therefore resort to an alterna- 
tive approach, and determine the effective bending rigid- 
ity by numerical integration of the mechanical equilib- 
rium equations in real space, dH/dcj) = 0. Specifically, 
we determine the end-deflection y{L) of a bundle of 4 fil- 
aments under a point force F at the distal end (s = L), 
given clamped boundary conditions at the proximal end 
(s = 0). The effective bending rigidity is then obtained 
from the expression, k^s = FL^/12y{L), and displayed 
in Fig. pU|) as a function of the crosslink shear stiffness, 
fcx, and a series of values for the pretwist, ujq- For sim- 
plicity we have assumed the filaments to be inextensible, 
ks ^ CO, and thus restricted ourselves to the decoupled 
and the intermediate regimes. 

For increasing pretwist the apparent stiffness decreases 
and asymptotically approaches the value without shear- 
stiffness. Thus, in pretwisted, helical bundles the 
crosslinks only have a limited ability to mechanically 
couple the different filaments together. The higher the 
twist the more the filaments act as if they were inde- 
pendent [2^ . The reason for this behavior is that in 
pretwisted bundles the filaments exchange their place and 
those that start on the top of the bundle soon are on the 
bottom. Thus, crosslink sites that would stay behind and 
lead to large shear displacements in untwisted bundles, 
can now catch up to make the effective shear deformation 
smaller. 



VI. CONCLUSIONS AND OUTLOOK 

We have presented a detailed study of the elastic and 
dynamic properties of bundles of semiflexible filaments 
(wormlike bundle model, WLB). It is found that a com- 
petition between the elastic properties of the filaments 
and those of the crosslinks leads to renormalized effective 
bend and twist rigidities that can become mode-number 
dependent. The strength and character of this depen- 
dence varies with bundle architecture, such as the ar- 
rangement of filaments in the cross section and pretwist. 

Two paradigmatic cases of bundle architecture have 
been discussed (see Fig. ^ : the first assumes filaments to 
be arranged homogeneously throughout the cross-section, 
for example on a square or triangular lattice. This ge- 
ometry is particularly relevant for F-actin bundles. The 
second architecture has the filaments arranged on the 
surface of a cylinder as is the case for microtubules. For 
all bundle architectures, the bending rigidity depends on 
mode-number qn as Ksin) kxq~'^. This is the shear 
dominated regime, as the bending rigidity is proportional 
to the shear stiffness, fcx, of the crosslinks [ll|]. It is 
in this parameter regime that the bundle behaves qual- 
itatively different than either a homogeneous beam (ob- 
tained in the "fully coupled" limit) or an assembly of 
"decoupled" filaments. Each architecture has its own 
universal ratio of maximal to minimal bending rigidity, 
independent of the specific type of crosslink induced fil- 
ament coupling. For microtubules (without pretwist) we 
find the ratio r = 1 + 2/ sin^(7r/A^), which is in reasonable 
agreement with the available experimental data. 

An important factor in determining the strength of 
crosslink-induced filament coupling is the pretwist (he- 
licity) of the bundle. Numerical computation shows that 
the effective bending rigidity decreases with increasing 
the pretwist. This has interesting consequences for mi- 
crotubules, where the amount of pretwist depends on 
the number of protofilaments, N. Different microtubule 
types are therefore predicted to have different variations 
in their effective bending rigidity. These predictions 
could be tested in experiments that are able to select 
the microtubule type and measure their bending rigidity 
independently. 

We have discussed several further observables, static 
and dynamic, that could be relevant to experiments. We 
have shown that the concept of the persistence length 
becomes ambiguous and depends on the observable used. 
The usual definition via the tangent-tangent correlation 
function is shown to lead to a persistence-length that 
depends on the scale of observation. Further observables 
that are eflFected by the mode-dependent bending rigidity 
are the force-extension relation or the buckling force [l2| . 
Interestingly, in the intermediate regime the latter is con- 
stant and independent of bundle length. 

The dynamic properties of bundles are characterized 
by a complex cross-over scenario which is in one-to-one 
correspondence with the three regimes of decoupled, in- 
termediate and fully-coupled bending. While decoupled 
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and fully- coupled bending display the usual t^/"* in the 
correlation function, it is shown that the shear-coupling 
leads to an intermediate asymptotic regime, where the 
correlations only grow as t^/^. The response functions 
for longitudinal and transverse forces also reflect these 
different regimes. In contrast to the WLC, they are not 
just proportional to each other but show distinct fre- 
quency dependences. These findings may be relevant for 
microrheological experiments, with imbedded bead par- 
ticles directly coupling to the transverse bundle fluctua- 
tions. 

Possibilities for future studies may be to look into the 
effects of filament fracture or lattice defects. The elas- 
tic energy represents a harmonic approximation which 
should be extended to include nonlinear effects. Espe- 
cially for microtubules one may expected nonlinear ef- 
fects to play an important role in bundle mechanics. For 
example, it may be important to consider that protofil- 
aments in their unstressed state are not straight but 
bend radially outwards. Aditional complications could 
also arise from the fact that some microtubules are not 
transvcrsly isotropic as we have assumed here, but have a 
" seam" , where protofilaments are offset relative to their 
neighbors. Failure modes under axial compressive forces 
have been discussed in a model for microtubules that 
starts from a transversly isotropic shell theory [sl]. It 
would be interesting to compare the results of a gener- 
alized WLB model - to include cross-sectional deforma- 
tions ~ with their results for the critical buckling forces. 
One mode of failure, the ovalization of the microtubule 
cross-section (Brazier effect) may for example be taken 
into account by adjusting the cross-sectional coordinates 
with the help of an "ovalization parameter" . This would 
then have to be determined together with the other de- 
grees of freedom from the equilibrium equations. 
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First, we have to show that there are no terms in the 
shear Hamiltonian, Eq. ([3]), that would couple the dif- 
ferent Euler angles d,(j),ip. To this end we use Eq. ^ 
in Eq. ([3]) and specialize to a square cross-section. The 
resulting shear Hamiltonian can then be written as 



shear 
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(Al) 



where yi = b ■ {i + \) and Zj = b ■ {j + i). Here, each 
filament is labeled by the pair of indices (i, j), which de- 
notes its location in the ith row and the jth column of 
the square cross-section. The only terms that couple the 
different Euler angles are b'^ip'(t)J2ij ^^'^ ^'^''P'^J^ij Vi- 
These identically vanish because of the symmetric ar- 
rangement of the filaments, Zj = and yi = 0. 

The remaining calculations are performed in Fourier- 
space. For the transformation we use cos-modes, which 
are appropriate for bundles with pinned boundary condi- 
tions. Writing the M-dependent part of the Hamiltonian 
in matrix form (iH = i '^ki "^kAkiui + uibi one needs 
the following formula valid for a Gaussian integral 

I Y\_dukcxp i-^'^UkAkiui -'^uibA (A2) 

= exp I ^ UkAkiUi - ^ uibi I , 
\ kl I J 



where u is obtained from 

d{l3H)/duk = 0. 



(A3) 



Having found the solution u, we can finally bring 
Eq. (K2\i in the form of Eq. ^ to read off the effective 
bend and twist rigidities. 
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APPENDIX A: EXACT SOLUTION FOR 
BENDING AND TWIST RIGIDITIES 

In this appendix the effective bending and twist rigidi- 
ties are calculated exactly from the WLB Hamiltonian, 
Eqs. H]) - (HI). To this end the u- variables have to be 
integrated over to define an effective WLC Hamiltonian. 



1. Bending 

Here we solve Eq. (jASp for the case of bending of the 
square bundle. It proves useful to introduce the dimen- 
sionless crosslink shear stiffness a = kxL'^ /kgS'^. We can 
then write Eq. (jA3p as 

{qL fuk - ad^uk = . (A4) 

We also defined the discrete second derivative d^Uk = 
Uk+i ~2uk + Uk-i- Note, that we are working in Fourier- 
space so all quantities should carry an additional sub- 
script relating to the mode- number n. As different modes 
don't mix, there is no harm in dropping it for the mo- 
ment. 

At the outer edges of the bundle k = —M, M — 1 we 
find 

(qL)'^ 

UM-l - UM-2 = UM-2 + b9 , (A5) 
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and 

(qL)^ 

u-M+i-u-M^ u-M + b9. (A6) 

a 

The Eq. (|A4p is solved by Uk = Am^ + Bm!t , where 
m± = i (a; ± ^x'^ — 4) and we have defined x = 2 + 
[qLY/a. The constants A^B are adjusted to fulfil the 
boundary conditions Eqs. (|A5[) and (jA6[) . This solution 
for the axial stretching variable Uk is plotted in Fig. [6l 

To obtain the approximated bending rigidity of 
Eq. (|10p one first has to insert the assumption = 
Aii • (fc + 1/2) into the Hamiltonian and then minimize 
with respect to the single variable Aw. 

This way one finds 




which has to be reinserted into the Hamiltonian to yield 

Eq. uni). 



2. Twist 

In the case of pure twisting the same analysis can be 
done to calculate the effective twist rigidity. The differ- 
ence to the bending case is that now the full 2d boundary 
value problem has to be solved. The equation determin- 
ing the axial stretching of the filament indexed by 
is 

{qLfu^j - aAu^j = 0, (A8) 

where the operator A — df + dj is the two-dimensional 
version of the discrete Laplacian encountered above. The 
finite difference operator di is defined as dfuij = Ui+ij 4- 
Ui-ij — 2uij. As for the case of pure bending, here, the 
twist variable ip enters only via the boundary terms. 

The classic theory to calculate the twist rigidity of 
beams has been given by Saint- Venant [13]. In this 
approach the axial displacements are found by solving 
Laplace's equation Au(y, z) = on the appropriate do- 
main of the cross-section. We see that Eq. (|A8p reduces 
to the Laplace equation in the limit a ^ kx ^ oo, which 
is the reason why in Fig. [5] the continuum limit is ap- 
proached with increasing shear stiffness at fixed bundle 
size. The remaining difference stemming from the dis- 
cretization into N filaments represents only a small effect. 
In Saint- Venant theory it is well known that for rectan- 
gular domains two types of solutions appear, depending 
on the aspect-ratio of the cross-section. We illustrate 
the different symmetry properties of these solutions in 
Fig.im 




FIG. 11: Illustration of the two types of solution obtained for 
the stretching Uij of filament Left; square cross-section 

of 40 * 40 filaments. Right: 30 ★ 40 filaments 

APPENDIX B: RELATION BETWEEN THE ^ 
OF FILAMENT k AND THE OF THE 
CENTRAL LINE 

This appendix gives some details on the description of 
the bundle kinematic degrees of freedom. The goal is to 
relate the generalized curvatures Qa,k of filament fc, to 
the rict of the bundle central line. That these need not 
necessarily coincide is best illustrated by discussing an 
example. If the central line is twisted but not bent, fls 
remains as the only non- vanishing component of the cur- 
vature vector. The filaments themselves, however, twist 
and bend as they trace out a helical path with radius R. 
Their bending energy is KfL{Rrt^)^ , of fourth order in 

As a second example consider the bending of the cen- 
tral line, r^i = i/R, where R is the radius of curvature. 
The filaments that lie at a distance b away from the cen- 
tral line naturally have a different radius of curvature, 
R±b, and thus a different bending energy. The correc- 
tion is again of higher order and scales as fe^fif . It turns 
out, that all effects like the two just mentioned only con- 
tribute to higher order. To lowest order we will find that 

We define the vector rfe(s) — ro(s) + R/£(s) to point 
to filament k at arclength position s. The central line 
of the bimdle is thereby given by ro(s). The position 
of each filament in the cross-section is parametrized by 
Ra;(s) = Afcdi(s) + Bkd2{s), where Ak and Bk are con- 
stants independent of arclength s and deformation of the 
bundle. This, in particular, implies that in the reference 
state the filaments are always straight and untwisted. We 
also need the derivative of R/j with respect to s, 

R;,^ = (Bfcfii - Akn2)t + QsRi (Bl) 

where we have used the Frenct-Seret equations and de- 
fined R^ = Bkdi{s) — Akd2{s). With this the tangent 
to filament k is given by 

tk = ^((1 + Bkfli - Akn2)t + rfsRfe ) , (B2) 

with an appropriate normalization factor, N. One finds, 
that the filament tangent is parallel to the central-line 
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only if the twist vanishes, f^a = 0. For finite bundle 
twist, the filament tangent is rotated around the vector 
Rfe relative to the central tangent t. The magnitude of 
the rotation is Q^Rk and depends on the distance Rk = 
|Rfe| = I'R-kl of the filament to the central line. 

In order to derive expressions for the remaining two 
unit vectors di.fc and d2.k let us assume, for the time be- 
ing, that no other rotation is allowed that may re-orient 
the local frame of filament k relative to the central frame. 
With this assumption filaments are not allowed to twist 
relative to their neighbors as this would correspond to a 
rotation around t^. In the bundles we consider, this in- 
ternal twist motion can safely be neglected as crosslinks 
provide for permanent rigid inter-filament connections. 
As explained below, these internal twist modes may nev- 
ertheless be important at the time of bundle formation. 

The local frame can then be related to the central 
frame by 



di t = di - Sfefissinat 



Oibk^sf 



(B3) 



d2,fe - 62 + Ak^s COS at + 0{bkn3 f (B4) 

One can now express the generalized curvatures of the 
filaments in terms of the fla of the central line. For sim- 
plicity we here consider the case of pure twisting where 
ill = ^2 = 0. One then finds 



^3.k = (ds^di^k) ■ d2,fc 



(B5) 



L 



)) 



Note, that the difference in length between the central 
line and the individual filaments should be taken into ac- 
count in the paramctrization in terms of arclcngth, how- 
ever, only contributes to higher orders. More formally, 
ds/dsk = 1 + 0{bn). 

In addition to twist, filaments also bend 



The twist energy should thus be written as 



(B13) 



twist 



where we defined the moments of the filament twist dis- 
tribution V""" = Efe i''k /N- 

As explained above, in crosslinkcd bundles this inter- 
nal twist can be assumed to be quenched at the time 
of bundle formation. In this case we need not treat the 
'i/'fc as dynamical variables for the discussion of bundle 
deformation. A finite -0' nevertheless gives the bundle 
a certain helicity and imposes pretwist, as discussed in 
Section IVl 



APPENDIX C: DERIVATION OF SHEAR 
DEFORMATION 

To calculate the geometric part A of the shear defor- 
mation, an expression for the arclength mismatch be- 
tween the two points on the filament pair is needed. The 
general expression is 



As 



-r;i 



(Cl) 



where the first contribution. As, derives from the possi- 
bility that the two points do not correspond to the same 
point, s, on the bundle central line. 

By using Eqs. (|Bip . (|B2p and expanding one finds 



A = As 



1 



(r: 



(C2) 



d3,fc 

k 



and 



^2,k = (ds^di^k) ■ d 



3,fc 



(B6) 
(B7) 

(B8) 
(B9) 



which reduces to the well known expression for the cur- 
vature of a helix if Jig = 0. 

Finally, let us comment on what would happen if we 
did allow the filaments to twist individually, and relative 
to their neighbors. This filament twist can be described 
by an additional rotation, ipk{s) of the local unit vectors, 
di,fe,d2,fe around the tangent tfc. 

The new vectors are then given by 

ei^k = cosV'fcdi^fc -I- sinV)fcd2,fc (BIO) 
62, fc = - sin'0fcdi,fc cos-0fcd2,fc , (Bll) 

and the bundle twist, Q^, can be calculated as before 

^l3,k = ^3 + ^'k ■ (B12) 



where, we defined b/fe = R; — R^ pointing from filament 
k to filament I. As in the formulation of Hq, only the 
leading order terms have been accounted for. The (bj^ -t)- 
term only contains bending deformations. This may be 
seen by setting t = e^; = const appropriate for a pure 
twisting of the central line. Then Aik ~ hik ■ = 
as the vector hik lies within the cross-section, that is 
perpendicular to the tangent. 

The last term corresponds to the arclength-difference 
acquired between two filaments at different distance to 
the central- line {Rk ^ Ri). It is clear that the filament 
farther out has to go a longer distance, so its crosslinking 
sites will stay back in comparison to that of its neighbor 
closer in the center of the bundle. However, twist also 
produces shear deformation between filaments which lie 
at equal distance to the central line (i?; = Rk). This 
is embodied in the extra term As, which we treat now 
(also see Fig. H]). To derive an expression for the shear 
displacement in this case, assume that the two filaments 
lie, separated by a distance bik = |b/fe|, on the surface of 
a cylinder of radius R. Twisting the cylinder makes the 
filaments wind around it, each taking an angle f^ai? to the 
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cylinder axis. The shear deformation then simply is As = 
fl^Rbik- For arbitrary orientation of the filament pair wc 
have to write As — ti^-hik-, where tik is the tangent to the 
midline between the filament pair. In agreement to what 



has been said above, this mechanism does not contribute 
to the shear displacement when the filament pair (I, k) is 
connected by crosslinks in radial direction. In this case, 
tik J- hik, and the shear deformation As = 0. 
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